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Abstract: The channeling of the ion recoiling after a collision with a WIMP changes the 
ionization signal in direct detection experiments, producing a larger signal than otherwise 
expected. We give estimates of the fraction of channeled recoiling ions in Si and Ge crystals 
using analytic models produced since the 1960's and 70's to describe channeling and blocking 
effects. We used data obtained to avoid channeling in the implantation of dopants in Si 
crystals to test our models. 
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1. Introduction 

Channeling and blocking effects in crystals refer to the orientation dependence of charged 
ion penetration in crystals. In the "channeling effect" ions incident upon a crystal along 
symmetry axes and planes suffer a series of small-angle scatterings that maintain them in 
the open "channels" in between the rows or planes of lattice atoms and thus penetrate much 
further into the crystal than in other directions. Channeled incident ions do not get close to 
lattice sites, where they would be deflected at large angles. 

The "blocking effect" consists in a reduction of the flux of ions originating in lattice sites 
along symmetry axes and planes, due to large-angle scattering with the atoms immediately in 
front of the originating lattice site, creating what is called a "blocking dip" in the flux of ions 
exiting from a thin enough crystal as function of the exit angle with respect to a particular 
symmetry axis or plane. 
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Channeling and blocking effects in crystals are related because the non-channeled incident 
ions are those which suffer a close-encounter process with an atomic nucleus in the crystal, 
namely those which pass sufficiently close to a lattice nucleus to be deflected at a large angle. 
After a close-encounter collision the deflected ion acts as if it was "emitted" from a lattice site. 
Channeling is many times observed as a lack of ions (incident at a small angle tp with respect 
to a particular symmetry axis or plane) deflected at a large-angle which form a "channeling 
dip" in the outgoing flux as function of the incident beam angle tp. As pointed out first by 
Lindhard pi, when no slowing-down processes are involved the "channeling" and "blocking" 
dips should be identical, when compared for the same particles, energies, crystals and crystal 
directions. 

Channeled ions loose their energy to electrons. They penetrate distances much larger than 
the characteristic separation of atoms along the channels, thus they interact with hundreds 
or thousands of lattice atoms. For energies in the keV range and above, channeled ions 
penetrate distances of at least several 10's of nm (see Appendix A, where we use the Lindhard- 
Scharff |2[ |3| model of electronic energy loss to calculate the penetration length of ions) . These 
are distances much longer than the separation of atoms along the channels, which are similar 
to the lattice constant, i.e. approximately 0.5 nm for Si and Ge (see Appendix B). 

The potential importance of the channeling effect for direct dark matter detection was 
first pointed out for Nal (Tl) by Drobyshevski Q and by the DAMA collaboration Q. The 
prospect of a daily modulation of the dark matter signal in direct detection due to channeling 
was recently raised by Avignone, Creswick and Nussinov || in Nal. 

In this paper we compute the channeling fraction of recoiling ions in Si and Ge crystals 
as function of the recoil energy and temperature. Si and Ge crystals are used in several 
direct dark matter detection experiments, such as CDMS J?]], CoGeNT ||, Edelweiss |J, 
TEXONO 0], EURECA gT|, HDMS @ and IGEX (D§. In a companion paper 0] we 



introduced the general ideas and analytic models |l], |3|, [l5|, |T(], 17, 18, 19, 20, 21, ^2| that we 
use to describe these phenomena in the context of dark matter detection, and applied them 



to Nal (Tl). For the reader familiar with Ref. 14] we would like to clarify which are the main 



differences between the calculations in Ref. [14] and in the present paper, besides the crystal 



structure (see Appendix B). In this paper we use a different expression for the continuum 
potentials (see Eqs. 2.1 to 2.7), which leads to a different expression for the critical channeling 
distance for axial channels (see Eq. 2.15). We also use a different way of deriving the critical 
distance for planar channels (see Eqs. 2.18 to 2.23). 

2. Model of Channeling 
2.1 Continuum models 

There are different approaches to calculate the deflections of ions traveling in a crystal. In 



"binary collision models" the ion path is computed by a computer program (see Ref. 23] 
for one of the first ones) in terms of a succession of individual interactions, each with one 
of the atoms in the crystal. Crystal imperfections and lattice vibrations are thus easily and 
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correctly taken into account. In "continuum models", reasonable approximations are made 
which allow to replace the discrete series of binary collisions with atoms by a continuous 
interaction between a projectile and uniformly charged strings or planes. These models allow 
to replace the numerical calculations by an analytic description of channeling, and provide 
good quantitative predictions of the behavior of projectiles in the crystal in terms of simple 
physical quantities. This is the approach we use here. 

The analytical description of channeling phenomena was initially developed mostly by 
J. Lindhard [BJ and collaborators for ions of energy MeV and higher, and its use was later 
extended to lower energies, i.e. hundreds of eV and above, mostly to apply it to ion implan- 



tation in Si. For the low energy range, we found most useful the work of G. Hobler [22], 
who in 1995 and 1996 perfected and checked experimentally previous continuum model pre- 
dictions [EJ] for axial and planar channeling at energies in the keV to a few 100 keV range, 



developed to avoid channeling in the implantation of B, P and As atoms in Si crystals [25]. 
This approach must be complemented by determination of parameters through data fitting or 
simulations. Moreover, lattice vibrations are more difficult to include in continuum models. 
Since we use a continuum model, our results should in last instance be checked by using some 
of the many sophisticated simulation programs that implement the binary collision approach 
or mixed approaches (e.g. p6|). 

Our calculation is based on the classical analytic models developed in the 1960's and 
70's, in particular by Lindhard [|, |, [fi| 0, ||, ||, |lj ||] . The fact that the de Broglie 



wavelengths of ions in the keV energy range are of the order of ~ 0.01 pm (and smaller at 
higher energies), thus much shorter than the lattice constant of a crystal (~ 500 pm, see 
Appendix B), justifies using a classical treatment. We use the continuum string and plane 
model, in which the screened Thomas-Fermi potential is averaged over a direction parallel to 
a row or a plane. This averaged potential is considered to be uniformly smeared out along the 
row or plane of atoms, which is a good approximation if the propagating ion interacts with 
many lattice atoms in the row or plane by a correlated series of many consecutive glancing 
collisions with lattice atoms. We are going to consider just one row, which simplifies the 
calculations and is correct except at the lowest energies we consider, as we explain below. 
There are several good analytic approximations of the screened potential. Except when 



said otherwise, in this paper we use Moliere's approximation, following the work of Hobler 22] 



and Morgan and Van Vliet [17, pL8| . Moliere's approximations of the continuum potentials are 
more complicated and also somewhat better than Lindhard's expressions, which we used in 
our paper devoted to Nal Jl4[. Lindhard's expressions are easier to manipulate algebraically 
to obtain different quantities of interest. Still in this paper we used some expressions derived 
from Lindhard's form of the potentials. 

In Moliere's approximation the axial continuum potential, as a function of the trans- 
verse distance r to the string, is 

U(r) = {2Z 1 Z 2 e 2 /d) f(r/a) = E4, 2 J(r/a), (2.1) 

where E is the energy of the propagating particle and ipi is a dimensionless parameter defined 
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Z\ , Z<i are the atomic numbers of the recoiling and lattice nuclei respectively, d is the spacing 
between atoms in the row, a is the Thomas-Fermi screening distance, a = 0A685h(Z^ 2 + 
Z 2 1/2 )- 2 / 3 |3|, [HJ (1.225xl0~ 2 nm and 0.9296xl0" 2 nm for a Si ion in Si and a Ge ion in Ge 
respectively, see Appendix B) and E = Mv 2 /2 is the kinetic energy of the propagating ion. 
Moliere's screening function |l5| for the continuum potential is 

3 

/(£) = 5>^o(/^)- (2-3) 
i=l 

Here Kq is the zero-order modified Bessel function of the second kind, and the dimensionless 



coefficients cti and are on = {0.1,0.55,0.35} and = {6.0,1.2,0.3} [27|, for i = 1,2,3. 
The string of crystal atoms is at r = 0. 

In our case, E is the recoil energy imparted to the ion in a collision with a WIMP, 

|q| 2 , , 

E = i-^L, (2.4) 

and q is the recoil momentum. 

The continuum planar potential in Moliere's approximation fl5| , as a function of the 
distance x perpendicular to the plane, is 

U p {x) = (2vrnZ 1 Z 2 e 2 a) f p (x/a) = E^ 2 J p (x/a), (2.5) 

where the dimensionless parameter ip a is defined as 

2 27rnZiZ 2 e 2 a 
= ^ 1 ( 2 -6) 

and n = Nd pc h is the average number of atoms per unit area, where N is the atomic density 
and dpch is the width of the planar channel, i.e. the interplanar spacing (thus, the average 
distance of atoms within a plane is d p = 1/ y/Nd pc h) ■ The subscript p denotes "planar" and 

3 

i=i 

where the coefficients ol% and are the same as above. The plane is at x = 0. 

Examples of axial and planar continuum potentials for a Si ion propagating in a Si crystal 
and a Ge ion propagating in a Ge crystal are shown in Fig. 

The continuum model does not imply that the potential energy of an ion moving near 
an atomic row is well approximated by the continuum potential U. The actual potential 
consists of sharp peaks near the atoms and deep valleys in between. The continuum model 
says that the net deflection due to the succession of impulses from the peaks is identical to 
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Figure 1: Continuum axial (black) and planar (green/gray) potentials for (a) Si and (b) Ge ions, 
propagating in the <100> axial and {100} planar channels of a Si or Ge crystal respectively. The 
screening radii shown as vertical lines are asisi = 0.01225 nm and aceGe = 0.009296 nm (see Appendix 
B). 



the deflection due to a force —U'. This is only so if the ion never approaches so closely any 
individual atom that it suffers a large-angle collision. Lindhard proved that for a string of 
atoms this is so only if 

U"(r) < ±E, (2.8) 

where the double prime denotes the second derivative with respect to r. Replacing the 
inequality in Eq. by an equality defines an energy dependent critical distance r c such 
that r > r c for the continuum model to be valid. Morgan and Van Vliet ]l7j also derived a 
condition for axial channels, similar to Eq. but with the factor 8 replaced by 16. 

The breakdown of the continuum theory for a planar channel is more involved than for an 
axial channel because the atoms in the plane contributing to the scattering of the propagating 
ion are usually displaced laterally within the plane with respect to the ion's trajectory. Thus 
the moving ion does not encounter atoms at a fixed separation or at fixed impact parameter 
as is the case for a row. Morgan and Van Vliet 17] reduced the problem of scattering from a 
plane of atoms to the scattering of an equivalent row of atoms contained in a strip centered 
on the projection of the ion path on the plane of atoms. They then applied Eq. 2^ as the 
condition for planar channeling to the fictitious string defined in this way (more about this 
below) . 



2.2 The transverse energy 

Lindhard proved that for channeled particles the longitudinal component v cos eft of the ve- 
locity, i.e. the component along the direction of the row or plane of the velocity, may be 
treated as constant (if energy loss processes are neglected). Then, in the continuum model, 
the trajectory of the ions can be completely described in terms of the transverse direction, 
perpendicular to the row or plane considered. For small angle 4> between the ion's trajectory 
and the atomic row (or plane) in the direction perpendicular to the row (or plane), the so 
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called "transverse energy" 

E ± = E sin 2 4> + U ~ Ecp 2 + U (2.9) 



is conserved. In Eq. 2.S relativistic corrections are neglected. 

Let ri be the initial position at which the WIMP nucleus collision occurs, i.e. if rj > 
the recoiling nucleus was displaced with respect to its position of equilibrium in a crystal 
row when it collided with a WIMP. We call 4>i the angle of the initial recoil momentum with 
respect to the row of atoms and E the initial recoil energy of the propagating ion. Given 
these initial parameters, the issue of where to define E± arises. Namely, we define 

E ± = Esintf + U(r*), (2.10) 

but there are different possible choices for r* , the position at which to measure the potential 
U. In our case, the recoiling ion leaves an empty lattice site, thus it moves away from an 
empty lattice site in the potential generated by its neighboring lattice atoms. So the potential 
the recoiling ion moves through at the moment of collision is very small, and the recoiling ion 
conserves its momentum and direction of motion until it gets very near the nearest neighbor, 
a distance d away along the string. At this moment, it is at a distance 

r* = ri + d tan fa (2.11) 

from its nearest neighbor. Therefore, as we did in Ref. [HJ], we will make the approximation 



of defining the potential entering into Eq. 2.1C at this position r* . 



2.3 Minimum distances of approach and critical channeling angles 

The conservation of the transverse energy provides a definition of the minimum distance of 
approach to the string, r m i n (or to the plane of atoms x m \ n ), at which the trajectory of the 
ion makes a zero angle with the string (or plane), and also of the angle ip at which the ion 
exits from the string (or plane), i.e. far away from it where XJ ~ 0. In reality the furthest 
position from a string or plane of atoms is the middle of the channel, whose width we call 
d a ch for an axial channel (d pc h for a planar channel). Thus, for an axial channel 

E± = U(r min ) = Ei? + U(d ach /2). (2.12) 

We proceeded in two ways to define the axial channel radius (d ac h/2) for the axial channels we 
included in our calculation. We used the contour plots of the axial continuum potentials plot- 
ted in a plane perpendicular to the channels shown in Fig. 3 of the 1995 paper of Hobler |z| 
to read off the channel radius d a ch/2 of the <100>, <110> and <111> axial channels in 
terms of the lattice constant a; ai . They are 0.25 a/ ai , 0.375 a/ ai , and \/0.2 2 + 0.12 2 &\ at = 
0.233 a/ a t, respectively. For the other axial channels we considered, <211> and <311>, we 
define the channel width (f ac h in terms of the interatomic distance d in the corresponding row 
as dach = l/y/Nd, where is iV the atomic density. For a planar channel we replace the axial 



potential at the middle of the axial channel U(d ac h/2) in Eq. 2.12 by the planar potential 
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at the middle of the planar channel U p (d pc h/2) (the channel width d pc h was defined after 
Eq. P). 

For axial channeling Lindhard equates the condition for channeling with the condition 
in Eq. for the validity of the continuum model. Replacing the inequality in Eq. by 
an equality defines an energy dependent critical distance r c , so that channeling can happen 



only if the propagating ion always keeps a distance r > r c . Morgan and Van Vliet [IT] use 



5 instead of 8 in Eq. 2.5, because this agrees better with their simulations of channeling in 
copper crystals. Following Hobler p3], we use here Morgan and Van Vliet's equation to define 



i.e. 



U"(r c ) 



d 2 



E, 



(2.13) 



With Moliere's form of the potential it is not possible to solve analytically for r c . Morgan 
and Van Vliet [|l7]] gave the following approximate analytical solution for the axial channeling 
minimum distance of approach, 



(2/3)aV^ [1 - (Va/19) + (a/700)] 



(2.14) 



with a = (Z\Z2e 2 d/a 2 E). This solution is not correct at low energies (high values of a). As 




CI « 



200 




Figure 2: Comparison of the exact numerical solution (solid black) of Eq. 2.13 for the critical distance 



of approach r c (E) and the approximate analytic expression in Eq. 2.15 (dashed green) as a function 



of y/a = \J ZiZ2e 2 d/a 2 E for (a) the high y/a (low energy) range, and (b) the low \/a (high energy) 
range. The Morgan and Van Vliet approximation to r c (E) in Eq. 2.14| is also shown (solid red- labeled 
MV). 



can be seen in Fig. g (and also in Figs. 8 and 13 of the 1995 paper of Hobler j22j) the steep 
increase in the approximate Morgan and Van Vliet solution at low energies (see the curve 
labeled "MV" in Fig. ^.a) is not present in the numerical solution (see the curve labeled 
"Exact" in Fig. ^.a) of r c . Instead of Eq. |2.14| we use here a better approximate analytic 



solution obtained by fitting a degree nine polynomial to the exact solution of Eq. 2.12, 

r c = a [0.57305Va-0.0220301(Va) 2 + 0.000728889(\/a) 3 

-0.0000155189( v / a) 4 + 2.04162 x lO -7 ^) 5 - 1.65057 x 10~ 9 (V«) 6 

+7.9749 x 10~ 12 (Va) 7 - 2.11041 x 10" 14 (Va) 8 + 2.35121 x lO" 17 ^) 9 ]. (2.15) 
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Eq. |2.15| is valid from E of 1 keV to 29 TeV ( which corresponds to values of y/a between 180 
and 0.000158). Fig. § shows a comparison of the exact numerical solution r c (E) of Eq. |2.8| 



and the approximate analytic solution Eq. 2.15 as a function of sfa. (divided by the screening 
distance a). The high and low \fa range in Fig. ^.a and |2|b respectively corresponds to 
low and high energies. The maximum percentage error between the exact solution and the 
analytic approximation we use is 11.5 %. 

Fig. U shows the critical distance of approach r c (E) in Eq. [2.15 as a function of energy 
of the propagating ion for several axial channels, for Si ions propagating in a Si crystal and 
Ge ions propagating in a Ge crystal. 

Since r c is the smallest possible minimum distance of approach to the string of a channeled 



propagating ion for a given energy E, i.e. 
monotonically with increasing r, then 



> r. 



U(r min ) < U(r c ). 



and the potential U(r) decreases 



(2.16) 



Using Eq. 2.12, this can be further translated into an upper bound on E± and thus on ip, the 
angle the ion makes with the string far away from it, 



V> < ME) 



U(r c (E)) - U{d ach /2) 



E 



(2.17) 



ip c {E) is the critical channeling angle for the particular axial channel, i.e. it is the maximum 
angle the propagating ion can make with the string far away from it (in the middle of the 
channel) if the ion is channeled. 
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Figure 3: Critical channeling distance of approach r c (E) in Eq. |2.14 as a function of energy of the 
propagating ion for several axial channels, for (a) Si ions propagating in a Si crystal and (b) Ge ions 
propagating in a Ge crystal. 



The critical distance r c (E) increases as E decreases (see Figs. ||, [| and |7| to 13). At low 
enough E, r c (E) becomes close to the radius of the channel d acn /2, and the critical angle 
ipc(E) (which is the maximum angle for channeling in the middle of the channel) goes to 
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zero (see Figs. @ to |l0| and ||, |l|). This means that there is a minimum energy below which 
channeling cannot happen, even for ions moving initially in the middle of the channel. This is 
a reflection of the fact that the range of the interaction between ion and lattice atoms increases 
with decreasing energy and at some point there is no position in the crystal where the ion 
would not be deflected at large angles. The existence of a minimum energy for channeling 



was found by Rozhkov and Dyuldya [28| in 1984 and later by Hobler [22] in 1996. It is clear 
that to compute r c {E) when it is not small with respect to the radius of the channel d ac h/2, 
and thus to compute the actual minimum energy for channeling, we would need to consider 



the effect of more than one row or plane (as done in Refs. [28] and [22]), thus our results are 
approximate in this case. 

For planar channeling we will follow the procedure of defining a "fictitious row" introduced 
by Morgan and Van Vliet [|l7], [22]| . They reduced the problem of scattering from a plane of 
atoms to the scattering of an equivalent row of atoms contained in a strip of width 2R (R is 
defined below) centered on the projection of the ion path onto the plane of atoms, and took 
the average area per atom in the plane, 1/Nd pc h to be 2R times the characteristic distance d 
between atoms along this fictitious row, i.e. 



d= l/(Nd pch 2R). 



(2.18) 



Once the width 2R of the fictitious row is specified, one uses the channeling condition for the 



continuum string model, Eq. 2.8, with the average atomic composition of the plane. For R, 
Morgan and Van Vliet used the impact parameter in an ion-atom collision corresponding to 
a deflection angle of the order of "the break-through" angle y/U p (0)/E. This is the minimum 
angle at which an ion of energy E must approach the plane from far away (so that the initial 
potential can be neglected) to overcome the potential barrier at the center of the plane at 
x = (namely, so that E± = U p (0)). For small scattering angles, the deflection angle 5 is 
related to the impact parameter, in this case R, as (see e.g. Eq. 2.1' of Lindhard [|]]) 



2E5 = -d U'(R), 



(2.19) 



where U' is the derivative of the axial continuum potential, and Morgan and Van Vliet define 
R by taking 5 = \/U p (0) /E. Using the Moliere's approximation for the potentials, Morgan 
and Van Vliet found the following expressions for R 



R 



which lead to the d value 



MV 



In B ZiZ 2 e 2 /aJEU p (0) 



d 



MV 



A aNd pch ln B Z 1 Z 2 e z /a J EU p (0) 



(2.20) 



(2.21) 



with coefficients A = 1.2 and B = 4. Morgan and Van Vliet [17] found discrepancies with 
this theoretical formula in simulations of binary collisions of 20 keV protons in a copper 
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crystal and adjusted the coefficients to A = 3.6 and B = 2.5. Hobler |22[] used both sets of 
coefficients and compared them with simulations and data of B and P ions propagating in Si 
for energies of about 1 keV and above. Hobler concluded that the original theoretical formula 
was better in his case (although Hobler proposed yet another empirical relation to define d). 



While Eq. 2.19 seems to provide a good condition for R, there is a channel dependent energy 
upper limit of applicability of its approximate analytical solution in Eq. [2.21 , because the 
logarithm in d MV approaches zero as E approaches {^Z\Z2e 2 / a) 2 /U P {Q). Close to this value 
of E there is an unphysical fast increase in d MV (and consequently in x c (E)) that indicates 
the break-down of the approximate solution d MV in Eq. |2.2l| (and, as shown in Fig. 13 of the 



1995 paper of Hobler [22|, is not found in other expressions of x c ). 

We decided to keep the Morgan and Van Vliet definition for R in Eq. |2.19| and use the 
following approximate analytical solution obtained by fitting a degree five polynomial in In y 



to the exact numerical solution of Eq. 2. IS 



R = a (0.716014 + 0.510922 In y + 0.12047(ln yf + 0.0180492(ln yf 
+0.00442459(ln y) A - 0.000824744(ln yf) , 



(2.22) 



where y = ZiZ 2 e 2 /a- s /EU p (0). 

Fig. shows a comparison of the exact numerical solution of Eq. 2. IE for R and its 



analytical approximation in Eq. 2.22] (divided by a) as a function of y. Also the approximate 
expression of Morgan and Van Vliet in Eq. [2.21 is shown in Fig. |] (labeled MV). The high 
and low y ranges in Fig. |].a and b respectively corresponds to low and high energies. The 
approximate solution is not valid at y < 0.15 which corresponds to E > 50 MeV for Si, 
and E > 700 MeV for Ge. Within its range of validity, the percentage error of the analytic 



approximation in Eq. 2.22 is less than 



ft! I a 




250 




Figure 4: Comparison of the exact solution (solid black) of Eq. 2.19| for R/a and its analytical 
approximation in Eq. 2.22) (dashed green) as a function of y — Z^Zie 2 / ' a^J EU p (0) for the (a) high y 
(low E) range and the (b) low y (high E) range. Also the Morgan and Van Vliet approximation to 



R/a in Eq. 2.20 is shown (solid red- labeled MV). 



Let us call r c {E) the critical distance obtained from Eq. 2.15 for the fictitious row, whose 
interatomic distance is d in Eq. 2. 1£ in which the distance R is given in Eq. m Then > the 
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minimum distance of approach for planar channeling is 

x c (E) = f c (E). (2.23) 

Fig. [| shows the plot of x c {E) (obtained from using Eq. 2.22| for the fictitious string) as a 
function of energy for the most important planar channels, i.e. {100}, {110}, {111}, {210} 
and {310}. Fig. || shows that we can safely extend our approximation to 50 MeV for Si ions 
in a Si crystal and to 700 MeV for Ge ions in Ge crystals. 
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Figure 5: Critical channeling distances x c {E) in Eq. |2.23| as a function of the energy of the propagating 
ion for different planar channels, for (a) Si ions propagating in a Si crystal and (b) Ge ions in Ge. 



Writing equations equivalent to Eqs. |2.12| and 2.16| for planar channels, namely 

E± = U p (x min ) = E{^) 2 + U p (d pch /2) 



and 



U p (x min ) < U p (x c (E)), 



(2.24) 



(2.25) 



we obtain an equation similar to Eq. 2.17 but for the maximum planar channeling angle, the 
critical planar channeling angle 



U p (x c (E)) - U p (d pch /2) 
E 



(2.26) 



For very small energies, for which x c (E) > d pc h/2 no channeling is possible (the maximum 
distance to any plane cannot be larger than half the width of the channel separating them) 
and ipc=0 (see Figs. ||, |?] to 10 and 14. b, 15. b). When x c {E) approaches the middle of the 
channel the effect of other planes should be considered, so our approximation of using the 
potential of only one plane is not correct in this regime. 

The static lattice critical distances presented in Figs || and |5| (also in the left panels of 
Figs. |^, ||, H and 10) do not include thermal effects. These are important and must be taken 
into account. They increase the critical channeling distances and consequently decrease the 



critical channeling angles as the temperature increases (as clearly shown in Fig. 12 and |l 
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Figure 6: Temperature dependent Debye model one dimensional rms vibration amplitude u\{T) 



(Eq. 2.27) of the atoms in (a) a Si crystal and (b) a Ge crystal. For comparison the Thomas-Fermi 
screening distances for two Si atoms and two Ge atoms, asiSi and aceGe respectively are also indicated 
(see App. A). 



2.4 Temperature dependent critical distances and angles 

So far we have been considering static strings and planes, but the atoms in a crystal are actu- 
ally vibrating. We use here the Debye model to take into account the zero point energy and 
thermal vibrations of the atoms in a crystal. The one dimensional rms vibration amplitude 
u\ of the atoms in a crystal in this model is [j~5|, |2l|| 



m{T) = 12.1 A 



1/2 



(2.27) 



where the 1/4 term accounts for the zero point energy, M is the atomic mass in amu, and 
T are the Debye temperature and the temperature of the crystal in K, respectively, and &{x) 
is the Debye function, 

*/ n 1 r tdt 

The Debye temperatures of Ge and Si are respectively 9 = 290 °K and 9 = 490 ° K g|, ||. 
The vibration amplitude u\ as a function of the temperature T is plotted in Fig. ^ for Si and 
Ge crystals. At room temperature (20 °C), u\ = 0.00849 nm for Ge and u\ = 0.00827 nm for 
Si. 

In principle there are modifications to the continuum potentials due to thermal effects, 
but we are going to take into account thermal effects in the crystal through a modification of 
the critical distances found originally by Morgan and Van Vliet |17| and later by Hobler |2^] 
to provide good agreement with simulations and data. For axial channels it consists of taking 
the temperature corrected critical distance r c (T) to be, 



r c (T) = Vrm + ldMT)] 2 , (2.29) 
where the dimensionless factor c\ in different references is a number between 1 and 2 (see e.g. 



Eq. 2.32 of Ref. [18] and Eq. 4.13 of the 1971 Ref. flj). 
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Figure 7: (a) Static critical distances of approach and Debye one dimensional rms vibration amplitude 
til of the atoms in the crystal at 20 °C and (b) critical channeling angles at 20 °C with temperature 
effects computed assuming c\ — C2 = c and c = 1 or c = 2 as indicated, as a function of the energy of 
propagating Si ions in the <100> axial (black) and {100} planar (green/light gray) channels of a Si 
crystal. 
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Figure 8: Same as in Fig. (?] but for the <110> axial and {110} planar channels of a Si crystal. 



For planar channels the situation is more complicated, because some references give a 



linear and others a quadratic relation between x c (T) and u\. Following Hobler 22] we use an 
equation similar to that for axial channels, namely 



x c (T) = ^xl{E) + {c 2Ul {T)]\ 



(2.30) 



where again C2 is a number between 1 and 2 (for example Barret [23] finds C2 = 1.6 at high 
energies, and Hobler |^2| uses C2 = 2). We will mostly use c\ = C2 = 1 in the following, to try 
to produce upper bounds on the channeling fractions. 

Using the temperature corrected critical distances of approach r c (T) and x c (T) (Eqs. |2.29 
and 2.30| ) instead of the static lattice critical distances r c and x c (Eqs. |2.15 and |2.23[) , in the 
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Figure 9: Same as in Fig. but for Ge ions propagating in the <100> axial (black) and {100} planar 
(green/light gray) channels of a Ge crystal. 
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Figure 10: Same as in Fig. [|but f° r the <110> axial and {110} planar channels. 



definition of critical angles, Eqs. [2.17 and 2.2C, we obtain the temperature corrected critical 
axial and planar angles, examples of which are shown in the right panels of Figs. ^ to 1C 
(ci = C2 = c and c = 1 or c = 2 at room temperature) . 

As shown in Fig. 11, with this formalism and using c\ = C2 = 2 we fit relatively well 
the critical angles measured at room temperature for B and P ions in a Si crystal (shown in 
green, or gray if color not available) in several channels, for energies between 20 keV and 600 
keV that Hobler [221 extracted from thermal wave measurements. 



Figs. 12 and 13 show clearly the temperature effects in the critical distances and angles 
for a specific channel, the <100> axial channel of a Si crystal and for a propagating Si ion. 
At small energies the static critical distance of approach is much larger than the vibration 
amplitude, so temperature corrections are not important. For small enough energies the 
critical distance becomes larger than the radius of the channel indicating that nowhere in the 
channel an ion can be far enough from the row of lattice atoms for channeling to take place 
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Figure 11: Comparison of theoretical (black lines) temperature corrected critical angles (with 
c\ = C2 = 2) and measured critical angles at room temperature extracted from thermal wave mea- 
surements (22j (green, or gray if color not available, dots joined by straight lines to guide the eye) as 
a function of the energy of (a) B ions and (b) P ions propagating in a Si crystal at T=20 °C, for the 
indicated axial and planar channels. 
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Figure 12: Static (green) and temperature corrected with c\ = C2 = c = 1 (black) (a) critical 
distances of approach (and u\(T) in red) and (b) the corresponding critical channeling angles as a 
function of the energy of propagating Si ions in <100> axial channels of a Si crystal. 
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(thus the critical channeling angle is zero). The exact calculation of the energy at which this 
happens would require considering the effect of more than a single row of atoms (which we 
do not do here) thus our results at these low energies are only approximate. As the energy 
increases, the static critical distance of approach decreases and when it becomes small with 
respect to the vibration amplitude u%, the temperature corrected critical distance becomes 
equal to (c\U\) which is larger for larger values of c\. When u\{T) becomes important in 
determining the critical distance, this becomes larger, and therefore the critical channeling 
angle become smaller, for higher temperatures. 

Figs. 14 and 15 show how the critical channeling angles change with temperature for four 
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Figure 13: Same as Fig. [l^ but using c\ — C2 — c — 2 in the temperature corrected critical distances 
of approach. 



particular channels, the <110> and <100> axial and the {110} and {100} planar channels, 
for Si ions in Si and Ge ions in Ge, respectively. In both cases the axial channeling angles are 
larger than the planar critical angles. The <110> and {111} critical channeling angles are 
the largest among the axial and planar channels respectively. For example, at E = 200 keV 
for Si ions in Si, the channels with the largest channeling angles are (in order of decreasing 
channeling angles): <110>, <100>, <211>, <111>, {111}, <311>, {110}, {100}, {310}, 
and {210}. We can clearly see that the critical angles become zero at low enough energies 




Figure 14: Temperature corrected critical channeling angles for T=40 mK, T=20 °C, T=600 °C, 
and T=900 °C as a function of the energy of a Si ion propagating in the (a) <110> axial channels 
and (b) {110} planar channels of Si crystal. 



(for which the critical distance of approach needed for channeling should be larger than the 
radius of the channel) indicating the range of energies for which no channeling is possible. 
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3. Channeling of recoiling lattice nuclei 

The channeling of ions in a crystal depends not only on the angle their initial trajectory makes 
with rows or planes, but also on their initial position. The nuclei recoiling after an interaction 
with a WIMP start initially from lattice sites (or very close to them), thus blocking effects 
are important. In fact, as argued originally by Lindhard in a perfect lattice and in the 
absence of energy-loss processes the probability of a particle starting from a lattice site to be 
channeled would be zero. This is what Lindhard called the "Rule of Reversibility. " However, 
any departure of the actual lattice from a perfect lattice, for example due to vibrations of 
the atoms in the lattice, violates the conditions of this argument and allows for some of 
the recoiling lattice nuclei to be channeled. Lattice vibrations are more important at hight 
temperatures and they have two opposite effects on channeling fractions: the probability of 
finding the atom which collides with a dark matter particle further out of its equilibrium 
position increases with increasing temperature thus channeling fractions increase, but the 
range of angles the trajectory of the propagating ion must make with the direction of the 
channel decreases with increasing temperature, which decreases the channeling fraction. 
We now estimate the channeling fraction using the formalism presented so far. 

3.1 Channeling fraction for each channel 

As in Ref. fj] we use a Gaussian function for the probability distribution g(r) for the per- 
pendicular distance r to the row at which the atom that collides with a WIMP is located at 
the moment of the collision due to thermal vibrations in the crystal 

g(r) = ^e W (-r 2 /2ul). (3.1) 



The one dimensional vibration amplitude u\ is given in Eq. 2.27. As explained in detail in 



Ref. [14], the channeled fraction Xaxial(-E') q) of nuclei with recoil energy E moving initially 
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in the direction q making an angle <j) with respect to an atomic row is given by the fraction 
of nuclei which can be found at a distance r larger than a minimum distance j-j m j n from the 
row at the moment of collision, determined by the critical distance of approach 

»oo 

„2 



drg{r) = exp (-r- , min /2u 1 ). 



It can easily be seen using Eqs. 2.10 , 2,11 , [2.12| and 2. IE that, because t7(rj + dtani 
if <f> > ipc no channeling can occur and Xaxi<d{E, 4 1 ) = 0- 
Using the condition 

E sin 2 <\> + Ufa + dtan 0) = U(r min ) < U(r c (E)), 

that implies the equality for the minimum initial distance rj m j n , 

U(r iMn + dttm(j)) = U{r c {E)) - E sin 2 0, 



(3.2) 

> U(r c ), 

(3.3) 
(3.4) 



in Ref. [14] we derived the following analytic expression for r, 
mation to the potential: 

Co 



1 + ^ 



exp 



-2 sin 2 (p/ipl) - 1 



?.,mm from Lindhard's approxi- 

(3.5) 



dtan ( 



where C is a constant, which was found experimentally to be C ~ v3 |Q]. We use here 
this equation because it is not possible to find a similar analytic expression using Moliere's 
approximation to the potential (although following Hobler we use Moliere's approximation 
to obtain the critical distances and angles). Tj im i n is a function of the temperature too, 
through r c (T). Notice that a small change in the critical distance r c {T) and thus in rj. m i n is 
exponentially magnified in the channeling fraction Xaxial (Eq. |3.2| ). This constitutes the most 
important difficulty to evaluate channeling fractions. The same happens for planar channels. 

For a planar channel, the Gaussian thermal distribution for the planar potential is one- 
dimensional (the relevant vibrations occurring perpendicularly to the plane), 

g{x) = (27ru 2 )~ 1/2 exp(-j; 2 /2Mf). (3.6) 

This is normalized to be 1 for — oo < x < +oo. In our calculations we only consider positive 
values of x, < x < +oo, for each plane, thus we multiply g(x) by a factor of 2 to find the 
fraction of channeled nuclei for a planar channel, 



Xplanar (E, 1 



2 g(x)dx = —= 



exp(— x 2 /2uf ) 



\/2' 



dx = erfc 



Ul 



X 'i i 



Ul 



where the minimum initial distance (derived in Ref. |l4j using Lindhard's planar potential) 
is 



(a/2) < 


C 2 - 


V(^ 2 /a 2 ) + C 2 - (x c /a) - (sin 2 0/^ 2 ) 


1 




V^/a 2 ) + C 2 - (x c /a) - (sin 2 0/^ 2 ) 





dp tan <^. (3.7) 
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Here <fi is the angle q makes with the plane, defined as the complementary angle to the angle 
between q and the normal to the plane, or as the smallest angle between q and vectors lying 
on the plane. Also in this case, Xpianarl-E, 4>) = if 4> > ipc- Note that Xi^ m \ n is also a function 
of the temperature through its dependence on x c (T). 



Using either Eq. |3.2| or Eq. 3.7, for an axial and a planar channel respectively, we define 
channeling fraction Xk f° r each channel k, which depends on the initial energy E, initial angle 
<j) and temperature T. Then we will sum over all channels and angles to obtain the total 
channeling fraction as function of E and T. 




Figure 16: Geometric channeling fraction Xicc{E, q) including 74 channels (Eq. 3.8) for each direction 



q plotted using the HEALPix pixelization of a sphere, for (a) a 200 keV Si ion recoil in a Si crystal, 
and (b) a 1 MeV Ge ion recoil in a Ge crystal, at 20 °C. Temperature effects in the lattice were 
included with c\—ci = \. The light green, light blue, dark blue, pink, red, and yellow colors indicate 
a channeling fraction of 0.5, 0.013, 7.5 x 10~ 4 , 4 x 10~ 5 , 10~ 5 and zero, respectively. 



3.2 Total geometric channeling fraction 

The geometric channeling fraction is the fraction of recoiling ions that propagate in the 1st, 
or 2nd, or ... or 74th channel. Here "geometric" refers to assuming that the distribution 
of recoil directions is isotropic. In reality, in a dark matter direct detection experiment, the 
distribution of recoil directions is expected to be peaked in the direction of the average WIMP 
flow. Here we examine this geometrical channeling fraction, and postpone the case of a WIMP 



wind to another paper [3C]. 



We include in our calculation only the most important channels, the same considered by 
Hobler Q. These are the <100>, <110>, <111>, <211> and <311> axial channels and 
the {100}, {110}, {111}, {210} and {310} planar channels. These constitute a total of 74 
channels, as explained in Appendix B. 

The probability Xrec(E, q) that an ion with initial energy E is channeled in a given 
direction q is the probability that the recoiling ion enters any of the available channels, i.e. 

Xrcc(£, q) = P{A X oiA 2 or ... or A u ). (3.8) 
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We compute this probability in the same way we did inRef. @, using a recursion of the addi- 
tion rule in probability theory and treating channeling along different channels as independent 
(see in Appendix C that this is a good approximation). 

Fig. [l^ shows the channeling probability Xrec{E, q) in Eq. |3.8| for a 200 keV recoiling Si 
ion in a Si crystal and a 1 MeV Ge ion recoil in a Ge crystal, at 20 °C. Temperature effects 
were included with c\ = oi = 1. The probability is computed for each direction and plotted 
using the HEALPix pixelization of a sphere. The light green, light blue, dark blue, pink, red, 
and yellow colors indicate a channeling probability of 0.5, 0.013, 7.5 x 10 -4 , 4 x 10~ 5 , 10~ 5 
and zero, respectively. 

To obtain the geometric total channeling fraction, we average the channeling probability 
Xrec(E,c[) over the directions q, assuming an isotropic distribution of the initial recoiling 
directions q, 

P TCC (E) = -L f X r CC (E,q)dn q . (3.9) 



4-7T _ 

This integral is computed using HEALPix [^] (see Appendix B of Ref. [14] for a complete 
explanation) . 

Our results for the geometric total channeling fraction for Si ions in a Si crystal and Ge 



ions in a Ge crystal are shown in Figs. 17, 18 and 19 for three different assumptions for the 



effect of thermal vibrations in the lattice, which depend on the values of the parameters c\ 
and C2 used in the temperature corrected critical distances of approach r c (T) and x c (T) in 
Eqs. 2.29] and 2.3C . The unrealistic case of assuming no vibrations in the lattice (except for 



vibrations of the colliding atom) corresponds to taking c% = C2 = and is shown in Fig [l7] 
for different temperatures because it provides an upper limit to the channeling fractions. In 
this case the channeling fractions reach a few % and they increase with temperature. 

In the literature, in other materials or for other channeling ions, values of c\ and C2 



between 1 and 2 are used. Thus, we show the c\ = C2 = 1 choice in Fig. 18 and the 



ci = C2 = 2 in Fig. 19. As the values of c\ and C2 increase, also the minimal distances from 
row or planes at which propagating ions must be to be channeled increase, thus the critical 
channeling angles decrease, what makes the channeling fractions smaller. If the values of c\ 



and C2 found by Hobler |22j and by us (see Fig. 11) to fit measured channeling angles for B 



and P ions propagating in Si apply also to the propagation of Si ions in Si, then the case of 



c\ = C2 = 2 in Fig. |19| should be chosen and the channeling fractions would never be larger 
than 0.3%. With c\ = c-i = 1 the channeling fractions reach about 1% and they increase with 
temperature. 

Please note that we have not considered the possibility of dechanneling of initially chan- 
neled ions due to imperfections in the crystal. Any mechanism of dechanneling will decrease 
the fractions obtained here. 

4. Main results and conclusions 

We have studied the channeling of ions recoiling after a collision with dark matter particles 
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Figure 17: Channeling fractions of (a) Si and (b) Ge recoils in a Si and a Ge crystal respectively, as 
a function of the ion energy for temperatures T=900 °C (orange or medium gray), 600 °C (green or 
light gray), 293 K (black), and 44 mK (blue or dark gray) in the approximation of c\ — C2 = ("static 
lattice"). This is an upper bound with respect to any non-zero values of ci and C2. Temperature effect 
are included in the vibrations of the colliding atom. 
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Figure 18: Same as Fig. [l?] but with c\ = ci = 1. 



within Si and Ge crystals. The calculations are similar because both crystals have the same 
structure. Channeled ions move within the crystal along symmetry axes and planes and suffer 
a series of small-angle scatterings that maintain them in the open "channels" in between the 
rows or planes of lattice atoms and thus penetrate much further into the crystal than in other 
directions. In order for the scattering to happen at small enough angles, the propagating 
ion must not approach a row or plane closer than a critical distance r c or x c respectively. 
These are given in Eqs. 2.15| and 2.23 for a "static lattice" (i.e. a perfect lattice in which 
all vibrations are neglected) and by Eqs. 2.2£ and 2.30 once temperature vibrations of the 
crystal lattice are taken into account. The temperature corrected minimum distances of 
approach (in Eqs. [2.29 and 2.30| ) depend on the one dimensional rms vibration amplitude 
ui(T) (Eq. |2.27| ), which increases with the temperature, through the coefficients c\ and C2- 
These dimensionless coefficients are found in the literature (for different ions and/or crystals) 
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Figure 19: Same as Fig. [I?] but with c\ = C2 = 2. 



to take values between 1 and 2. Channeled ions must have trajectories that at large distances 
from the atomic rows or planes must make an angle with respect to these rows or planes smaller 
than a critical angle given in Eqs. 2.17 and 2.26 for axial and planar channels respectively. 

The critical angles depend on the temperature through the minimum distances of ap- 
proach: as these increase with increasing temperatures, the critical angles decrease, what 
makes the channeling fraction smaller. However, there is a second temperature effect which 
makes the channeling fractions larger as the temperature increases: the vibrations of the 
atom which collides with the dark matter particle. Thus, the channeling fraction of recoiling 
ions is strongly temperature dependent. Depending on which of the two competing effects 
is dominant, the channeling fraction may either increase or decrease as the temperature in- 
creases. Increasing the temperature of a crystal usually increases the fraction of channeled 
recoiling ions (see Figs. 17 and|l^), but when the values of c\ and C2 are large (i.e. close to 
2) so the critical distances increase rapidly with the temperature, the opposite may happen 
(see Fig. ||). 

The vibrations of the atom colliding with the dark matter particle are essential to have a 
non-zero probability of channeling of the recoiling ion. A nucleus ejected from its lattice site 
by a collision with a dark matter particle is initially part of a row or plane. Thus, the recoiling 
nuclei start initially from lattice sites or very close to them. This means that blocking effects 
are important. In fact, as argued originally by Lindhard |]l|, in a perfect lattice and in the 
absence of energy-loss processes the probability of a particle starting from a lattice site to be 
channeled would be zero. This is what Lindhard called the "Rule of Reversibility. " However, 
vibrations of the atoms in the lattice violate the conditions of this argument and allow for 
some of the recoiling lattice nuclei to be channeled. The channeling fraction Xaxiab 

Eq. ETa 

or Xpianan Eq. 3.7 for axial and planar channels respectively, is given by the fraction of nuclei 
which can be found further than a minimum distance 

^2,min or Xi mm away from a row or plane 
at the moment of collision. This fraction increases as u\{T) increases. This is the effect that 
dominates the temperature dependence in Figs. 17 and Qq, in which the geometric channeling 
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fractions increase with increasing temperature. 

However, rj m ; n , Eq. BJjrf O-T 3^i,min; Eq. |3.7| , increase with increasing critical distances and 
this decreases the channeling fraction. The increase of the critical distances with temperature 
is more accentuated for large values of c\ and 02- This can be seen in Fig. JO], in which 
ci = C2 = 2 and some channeling fractions are larger at lower temperatures. 

The unrealistic case of assuming no vibrations in the lattice (except for vibrations of the 
colliding atoms) corresponds to taking c\ = C2 = (this is what we call the "static lattice" 
approximation) shown in Fig. 17 for different temperatures, provides an upper limit to the 
channeling fractions. This is the limiting case of the possibility that ci and C2 are smaller than 
1, in which the channeling fractions reach a few % at energies of 100's of keV and increase 
with temperature. 

We show the c\ = C2 = 1 choice in Fig. [l^ and the c\ = c<i = 2 in Fig. 19. If the values 
found by Hobler [^] and by us (see Fig. |TT| ) to fit the measured channeling angles for B and 
P ions propagating in a Si crystal apply also to the propagation of Si ions in Si, then the 
case of ci = C2 = 2 should be chosen and the channeling fractions would never be larger than 
a few 0.1%. In this case, as mentioned above, the channeling fractions at some energies are 
higher at lower temperatures. The c\ = c% = 1 case, instead, leads to maximum channeling 
fractions of roughly 1 %, which increase with increasing temperature. 

Notice that a small change in the critical distances r c (T) or x c (T) and thus in the initial 
minimum distances of approach rj jm ; n or Xj )m i n is exponentially magnified in the channeling 
fractions Xaxiab Eq. [T^, or Xplanar> Eq. [O]. This constitutes the most important difficulty to 
evaluate channeling fractions in the models we use. 

Notice too that we have not considered any mechanism of dechanneling of the channeled 
ions (due to irregularities in the crystals, for example) which would decrease the channeling 
fractions. 
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A. Penetration length of channeled ions 

Fig. 20 shows the maximum distance, x majX (E) a channeled ion with initial energy E propa- 
gates in a crystal channel, according to the Lindhard-Scharff [Q, [3| model of electronic energy 
loss, for a Si ion channeled in a Si crystal and a Ge ion in a Ge crystal. This model is valid for 
small enough energies, E < (Mi/2)Z^ 3 Vq (where vq = e 2 /h = 2.2 x 10 s cm/sec is the Bohr's 
velocity 0. Mi and Z\ are the mass and charge of the propagating ion) which is E < 24.3 
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MeV for a Si ion propagating in a Si crystal and E < 188.7 MeV for a Ge ion propagating in 
a Ge crystal. In this model the energy E(x) as a function of the propagated distance x and 
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Figure 20: Maximum distance x max (E) traveled by channeled Si ions in Si (black) or Ge ions in Gc 
(green/gray). 



the initial energy E is the solution of the following energy loss equation 

dE u 
-— = Kv, 

dx 



where v = w2EjM± is the ion velocity and K is the function 

£ e 8vre 2 Na Z 1 Z 2 



K 



zi+zi 



3 

2 \ 2 



(A.l) 



(A.2) 



Here £ e is a dimensionless constant of the order of Zf || , AT is the number of atomic centers 
per unit volume, ao — 0.53 A is the Bohr radius of the hydrogen atom. Explicitly, an ion 
with initial energy E at x = has energy 

E{x) = E ( 1 - (A.3) 



(A.4) 



after traveling a distance x. The range of the propagating ion is 

V^MiEn 



%max(E) 



K 



Fig. 20 shows that even at energies of a few keV a channeled ion interacts with hundreds 



of lattice atoms. The characteristic interdistance of atoms along the channels is the lattice 
constant, i.e. approximately 0.5 nm for Si and Ge crystals (see Appendix B). 
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B. Crystal structure of Si and Ge 



Silicon (Si) and Germanium (Ge) crystals have a diamond cubic type lattice structure which 
consists of two interpenetrating face centered cubic (f.c.c.) lattices, displaced along the body 
diagonal of the cubic cell by one quarter of the length of the diagonal. The unit cell, shown 
in Fig. 21, has 8 atoms. The lattice constant, the side of the cube in Fig. 21, is ai at = 0.5431 
nm for Si and 0.5657 nm for Ge (from the Table 3.4 of Ref. plf). 




Figure 21: Unit cell of a Si or Ge crystal (a) in three dimensions, and (b) projected on a plane. The 
black spheres represent Si or Ge atoms. 

The atomic mass and atomic number of Si and Ge are M$i = 28.09 amu, Mc e = 72.59 
amu, Zgi = 14 and Zq € = 32. The Thomas-Fermi screening distances for two Si atoms 
and two Ge atoms are asiSi = 0.4685A(Z^ 2 + Z«^ 2 ) -2 / 3 = 0.01225 nm and aceGe = 
0.4685A(Zg / e 2 + Z^)~ 2 / 3 = 0.009296 nm respectively. Once an origin of the coordinate 
system is fixed on a lattice point O, any position vector of a point on the crystal lattice can 
be written as R = nia + 722b + 713c with m, 712, and 713 specific integer numbers. The vectors 
a, b, and c are the basis vectors of the crystal lattice, and are three noncoplanar vectors 
joining the lattice point O to its near neighbors. For Si and Ge, the three vectors a, b, c form 
a Cartesian frame and their length is aiat/4. The integers m, 772, and 77.3 can be positive, 
negative, or zero. The direction of a crystal axis pointing in the direction R is specified by 
the triplet [71171272.3] written in square brackets, when 774, n%, and 773 are positive or zero. 
Note that if there is a common factor in the numbers n\, 772, 773, this factor is removed. 
Moreover, negative integers are denoted with a bar over the number, e.g. —1 is denoted as 1 
and the — y axis is [010] direction. The plane perpendicular to the [77177-2773] axis is denoted 
by (771772773). For example, the plane perpendicular to the [100] axis is denoted by (100), and 
that perpendicular to [101] by (101). The integers m, 772, and 773 are called Miller indices. 

In a cubic crystal, because of the symmetry of the unit cell, the directions [100], [010], 
[001], [100], [010] and [001] are equivalent. All directions equivalent to the [771772773] direction 
are denoted by <77i772T73> in angular brackets. For example, <100> indicates all the six 
directions mentioned. Similarly, <211> and <311> indicate twelve different directions each. 

When the unit cell has cubic symmetry, we can indicate all planes that are equivalent 
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to the plane (hkl) by curly brackets {hkl}. For example, the indices {100} refer to the six 
planes (100), (010), (001), (100), (010), and (001). The negative sign over a number denotes 
that the plane cuts the axis on the negative side of the origin. Similarly, <210> and <310> 
each indicate twelve different planes. 

Counting all the axes and planes we mentioned above, the total is 74. 

The interatomic spacing d in atomic rows and the interplanar spacing d pc h ("pch" stands 
for "planar channel") of atomic planes of monatomic diamond crystals, are obtained by mul- 



tiplying the respective lattice constant by the following coefficients [15]: 

• Rows: < 100 >: 1, < 110 >: l/y/2, < 111 >: 3\/3/4, < 211 >: y/6/2, < 311 >: 3\/lI/4 

• Planes: {100} : 1/4, {110} : l/2\/2, {111} : y/3/4, {210} : l/(4\/5), {310} : l/(2\/T0) 
The Debye temperatures for Si and Ge are = 490 K and = 290 K, respectively 

MM. 



C. Probability of correlated channels 



In this paper, as we did in Ref. [14], we treat channeling along different channels as inde- 



pendent events when computing the probability x rec (.E, q) i n Eq.3^ that an ion with initial 
energy E and direction q enters any of the available channels. Available channels are those 
whose axis or plane, respectively, form and angle with the direction q smaller than the critical 
channeling angle for the particular channel. 

In Appendix D of Ref. [jl4| we showed that we can obtain an upper limit to the channeling 
probability of overlapping channels by replacing the intersection of the complements of the 
integration regions in Eqs. [T^ and |3.7| with the inscribed cylinder of radius tmin equal to 
the minimum of the rj im i n or a^mm among the overlapping channels. We find that the two 



methods give practically indistinguishable results for Si and Ge (as we did in Ref. [14] for 
Nal), as clearly shown in Fig. ^2] for some particular examples. Thus, the method we use is 
adequate for our purpose of providing upper bounds to the channeling fractions. 

Fig. |23| shows the channeling fractions of Si ions propagating in a Si crystal and Ge ions 
propagating in a Ge crystal for individual channels with c% = oi = 1 and T= 293 K. The 
black and green (or gray) lines correspond to single axial and planar channels respectively. 



Fig. 23 shows that at low E channeling is dominated by axial channels which do not 
overlap, so treating them as independent is strictly correct. However, at the transition energy 
of 1 to 10 MeV at which axial and planar channels are both equally important, and at 
higher energies at which planar channels dominate, the overlap of one axial and two or 
more planar channels, or the overlap of two or more planar channels among themselves, 
makes the channeling along them not necessarily uncorrelated. Still we find that considering 
channeling along different channels as independent is a good approximation if we are interested 
in providing upper bounds to the channeling fractions. 
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Figure 22: Maximum channeling fractions for c% = = c and c = (dashed green), c = 1 (dashed 
yellow) and c = 2 (dashed cyan) compared with the results for the same models of our method of 
Section 3.2 (solid black lines) or (a) Si ions propagating in a Si crystal and (b) Ge ions propagating 
in a Ge crystal at T = 293 K. Notice that the lines for the same model overlap. 
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Figure 23: Channeling fractions of (a) Si ions propagating in a Si crystal and (b) Ge ions propagating 
in a Ge crystal for single planar (green/gray lines) and axial (black lines) channels, as function of the 
recoil energy E, for T= 293 K and c\ = c 2 = 1. 
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